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ABSTRACT 

Context. We introduce a new technique, called the Sliding Two-Dimensional Fluctuation Spectrum, used for detecting and character- 
ising the temporal changes of drifting subpulses from radio pulsars. The method was tested using simulated data as well as archived 
observations made with the Westerbork Synthesis Radio Telescope (WSRT) at wavelengths of 92 and 21 cm. 

Aims. The drifting subpulse phenomenon is a well known property of radio pulsars. It has been studied extensively since their 
discovery, however the properties of the temporal behaviour of drifting subpulses are not fully explored. The drifting can also be 
non-coherent and the presence of effects like nulling or drift rate changing can mask the drifting behaviour of the pulsar. The S2DFS 
is a robust method for investigating this phenomenon and by introducing it we aim to expand our knowledge of the temporal drifting 
subpulse properties. 

Methods. Our new analysis method uses horizonally collapsed fluctuation spectra obtained with the Two-Dimensional Fluctuation 
Spectrum (2DFS) method. Stacking the collapsed spectra obtained in a 256 pulse window which slides by a pulse at a time produces a 
map of the collapsed fluctuation spectrum (S2DFS plot). By analysing the maps one can easily determine the presence of any temporal 
drift changes. 

Results. Simulated data showed that the technique can reveal the presence of any temporal changes in drift behaviour like mode 
changing or nulling. We have also analysed data of three pulsars, PSRs B0031-07, B1819-22 and B1944+17, which were selected 
based on the quality of the data and their known drift properties. All three sources are known to exhibit mode changes which could 
easily be seen in the S2DFS. 

Conclusions. The results from the analysis of the data sets used in this paper have shown that the S2DFS method is robust and 
complimentary to the 2DFS method in detecting and characterising the temporal changes in drifting subpulses from radio pulsars. 

Key words, methods: analytical - pulsars: general 

1. Introduction The workof lBacked(ll970allRll973h and lBacker et al.l ( ll975h 

systematised the knowledge and introduced a better fluctuation- 
spectral analysis technique by applying the Discrete Fourier 
In the years follow ing the serendipitous discovery of pulsars Transform (DFT) to regions narrower than the pulse width, that 
(lHewish et al.119681) . the results from analysis of long sequences ^^^^^ ^^^^^ ^"'^f ^ . ^"^^ .^^ich insight into the proper- 
of pulses revealed characteristic "secon d periodic pulsa tions" ^les of subpulse modulation within a pulse profile. This tech- 
seen in the plots of the subpulse intensity (iDrake & Craf til 19681) . "^^l^^' ^no^n as the Longitude-Resolved Fluctuation Spectrum 
This phenomenon, called "class 2 period" (denoted by P2, see (LRPS), became very wide y used in the analysis of the subpulse 
Fig. D, expressed itself as a constant change of the phase of ^^^^ properties of pulsars. Apart from discovering and character- 
successively appearing subpulses within a fixed longitude range. i^^^S values, the LRFS method could not resolve informa- 
In order to investigate thi s new phenomen on, nu merous stud- ^1°" subpulses drift (i.e. value). This insformation was 
ies were conducte d (e.g. E ^ce & Craft l^; ' TavloretalJ "bt^i"^'! applying the Fourier transform to a series of pulses, 
[19691 ISuttonetal. 1970: £olsi970). The authors ^^^^ ^^^^^ '^P'^'^f ^ visual mspection of hig h sig^ 
performed fluctuation-spectral analysis using Fourier techniques to-""!^^ I'^^i" ^^/N) observations. Edwards & Stappers ( 2002D 
which revealed another type of characteristic frequency at which Presented a new technique, the Two-Dimensional Fluctuation 
a significant amount of power was present in the fluctuation Spectrum (2DFS) which could easily determine the P2 value, 
spectra of pulsars and provided evidence of a periodic modu- The 2DFS is an extension of the LRFS, in a sense, where, 
lation of pulse intensities. These intensity modulations were de- ^PP^^''^ ^FT along vertical lines of constant 
scribed by a pattern period P3 (see Fig. [D, which is the time P^^^^ in the pulse stack (the result of which gives the LRFS), 
between the successive appearance of subpulses at a particular ^ ^^'^0"^ ^FT is applied to each row of the complex version 
pulse phase, expressed in pulsar periods, P„ ( Sutton et al. 197(1). ^RFS, thus producing the two-dimensional power spec- 
Although very useful, these techniques had their limitations, and ^mm. While this results in the loss of pulse longitude informa- 
a better technique was needed. ^1°"' " P^vides the P2 value. The LRFS and 2DFS techmques 
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have proven to be robust and effective for discovering and char- 
acterising subpulse drifting from radio pulsars as shown in the 
recent work by Weltevrede et al. (2006, 2007). The Harmonic- 
Resolve d Fluctuation Spectrum techn ique, called HRFS, devel- 
oped by Deshpande & RankinI (1200 ll) . can also be used for sub- 
pulse drifting analysis and is mathematically equivalent to the 
2DFS. 

In their work, Weltevrede and collaborators examined a large 
sample of pulsars observed at two frequencies in order to study 
their subpulse modulation properties. A large fraction of all pul- 
sars (about 40% of a total sample of 185 sources) were re- 
ported to exhibit the drifting subpulse phenomenon. Those pul- 
sars whose P3 values were well defined have been classified 
as coherent drifters, while the sources for which P3 could not 
be well defined, were assigned to the diffuse class. Due to var- 
ious reasons like low S/N, interstellar scintillation effects, se- 
vere radio frequency interference during observations or condi- 
tions intrinsic to the observed source, not all pulsars in the sam- 
ple were observed to ex hibit the drifting subpulse phenomenon 
dWeltevrede et al.ll2006l) . 

Let us now consider the pulsar intrinsic effects: The first in- 
trinsic property that might prevent the detection of drifting is the 
viewing geometry. In this case the line of sight cuts the magnetic 
pole centrally and only longitude stationary subpulse modula- 
tion is expected. The second intrinsic property can be assigned 
to the refractive properties of the pulsar magnetosphere leading 
to a distorted pulsar signal. We will now focus on the other ef- 
fects potentially responsible for the lack of drifting: nulling and 
mode changing. 

Nulling characterises itself as a sudden diminishing or ces- 
sation of pulse energy below the level of dete ction for a cer- 
tain amount of time, and was first reported by iBackeij (Il970cl) 
in the observations of PSRs B0826-H06, Bl 133-1-16, B1237-H25 
and B1929-I-10. In Fig.[T]one can see an example of a null be- 
tween pulses 31 and 38 in the pulse stack of PSR B003 1-07 
Studies which followed this discovery (e.g . .Ritchingsl 119761: 
lRankinlll986l: lBiggslll992l: IWang et al.[|200% have shown that 
nulling can occur on random timescales while the duration can 
vary between one or more periods to longer timescales of hours 
to days. The "Nulling Fraction", which is the percentage of 
time for which the pulsar nulls ranges from almost zero (for the 
Vela pulsar ) up to 93 % (for PSR J1502-5653) as shown by 
I Wang et al] (|2007|) . It is also worth noting that although nulling 
is a broadband effect, nulls do not always occur si multaneously 
at different simultaneo usly observed frequencies (iDavies et al.l 
fT984tlBhat et"ani2007h . 

Pulsars are known to form extremely stable profiles made 
after averaging just a few hundred pulses dHelfand et al.lll9'75t 
iRathnasree & Rankin 1995). but in some cases they change 
the shape of their average profiles as first observed in PSR 
B1237-H25 and reported bv [Backed (Il970bl) . This behaviour, 
called mode changing, is known to happen i n about 12 pulsars, 
mostly with multicomponent ave rage profiles (iBartel et all 19821: 
iRankinI [19861: IWang et al]|2007h . The mode changing happens 
suddenly, it can include two or more modes and can last from a 
few to thousands of pulse periods. It is also important to note that 
mode changing is closely connected with the nulling and dr ifting 
subpulse phenomena. For example. iHuguenin et al.l (1 19701) have 
first shown that PSR B003 1 -07 changes its drift rate in each of 
its modes, while PSR B0809-I-74 after many or all nulls, emits i n 
a mode different from the normal one (van Lee uwen et al.l2002l) . 

The presence of the above mentioned effects can disrupt the 
subpulse drift of the pulsars and change their subpulse modu- 
lation periodicities which, depending on the magnitude of the 
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Fig. 1. An example of sequence of one hundred successive single 
pulses of PSR B003 1 -07 observed at 92 cm. The subpulses ap- 
pear earlier with increasing pulse number and are arranged into 
so-called drift bands. Two successive drift bands are vertically 
separated by f 3 and horizontally by ^2- 

effects and the S/N of the observation, can influence the fluc- 
tuation analysis. We introduce an extension of the 2DFS tech- 
nique called the S2DFS for discovering and characterising the 
temporal changes of drifting subpulses from radio pulsars. We 
begin with a more detailed description of the 2DFS technique 
for analysing subpulse modulation phenomena. We then intro- 
duce the new technique and apply it to the simulated and archival 
data of a few well known pulsars observed with the Westerbork 
Synthesis Radio Telescope (WSRT). Finally, we discuss the im- 
plications of the application of our method to pulsars with unsta- 
ble drifting behaviour, multiple drift modes or those affected by 
the nulling phenomenon. 

2. Methods of analysis 

2.1. 2DFS 

In order to analyse the subpulse modulation phenomena using 
the methods described in this work, for each data set, a so- 
called pulse stack has to be formed. To create a pulse stack, 
a de-dispersed time series is transformed into a longitude and 
time-resolved representation of the pulsar signal. This is done 
by stacking time samples from consecutive pulses, according to 
their phases, as calculated using an ephemeris and the TEMPO 
software packageQ. Fig. [Upresents an example of a pulse stack of 

' http://www.atnf.csiro.au/research/pulsar/tempo 
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one hundred pulses from PSR B003 1 -07 plotted one above the 
other, where the pulse longitude is plotted on the horizontal axis 
and the pulse number on the vertical axis. One can easily see 
that the subpulses appear earlier with increasing pulse number 
and are arranged into so-called drift bands. 

The basic method to estimate the presence of subpulse mod- 
ulation is to calculate the longitude -resolved modulation index 
(LRMI: lEdwards & Stappersll2002l) . The LRMI is the measure 
of the factor by which the intensity varies from pulse to pulse for 
each pulse longitude. In the right-hand plot in Fig.|2]the average 
profile and the LRMI are shown in the top panel for a simulated 
pulse stack. The solid line corresponds to the average profile, 
which is normalised to the peak intensity, while the solid line 
with error bars denotes the LRMI values. We only show LRMI 
values that are detected with a significance more than 3cr. 

Although the LRMI is the first evidence of subpulse mod- 
ulation, it does not indicate whether this modulation is system- 
atic or lacks organisation. To obtain such information the LRFS 
has to be calculated. This is done by dividing the pulse stack 
into blocks of 2" pulses and applying a DFT at each pulse lon- 
gi tude bin in each block (for deta ils of the analysis we refer 
to lEdwards & Stappersll2002l l2003h . The final spectrum is pro- 
duced by averaging each of the longitude-resolved spectrum in 
each block. If the pulsar exhibits subpulses which are periodi- 
cally modulated, then a region, so-called feature, of enhanced 
power will be visible as a dark region in the greyscale of the 
power in the LRFS. The vertical position of the feature is given 
in cycles per period (cpp) and corresponds to a period of mod- 
ulation, P3 expressed in pulsar periods Pq, while the horizontal 
position denotes the pulse longitude position at which the mod- 
ulation occurs. In the right-hand plots in Fig. |2] the LRFS of the 
simulated pulse stack is shown below the top panel. 

While the LRMI allows one to detect the presence of sub- 
pulse modulation, and the LRFS determines whether subpulse 
modulation is disorganised or periodic, the information about 
whether the subpulses are drifting and with what rate is not re- 
solved. Such information can be acquired by applying the 2DFS. 
The procedure is similar to the calculation of the LRFS where, 
except now after applying the DFT to the pulse stack along the 
vertical lines of constant phase, the DFT is applied again, but 
now across each row of the complex LRFS spectrum. In the 
right-hand plot in Fig. |2]the 2DFS for simulated pulse stack is 
plotted in the panel below the LRFS plot. The vertical axis of 
the resulting spectrum is the same as in the LRFS, but the hori- 
zontal axis is now given in cycles per period and corresponds to 
a horizontal separation of the drift bands (P2 expressed in pulse 
periods). If the feature is significantly offset from the vertical 
axis {P0/P2 + 0) it means that subpulses drift in a certain longi- 
tude range have a preferred drift direction. A negative or positive 
value of f 2 denotes the dkection of the drift, i.e. whether sub- 
pulses drift towards the leading or trailing edge of the pulsar's 
average profile respectively. The power of the 2DFS is vertically 
(between dashed lines) and horizontally integrated, resulting in 
the plots in the side and bottom panels, which are used for better 
understanding the structure of the feature. To measure values of 
f 2 and P3 one has to calculate the centroid of a rectangular re- 
gion in the 2DFS containing the feature (Eq. 6: IWeltevrede et akl 
[2OO6), while errors for those values are estimated by applying 
the same procedure to an equivalent rectangular region contain- 
ing only noise. 



2.2. S2DFS 

The 2DFS method is a very effective tool for subpulse drift 
detection and analysis. Its advantage over other methods lies 
in a process of averaging multiple fluctuation spectra, obtained 
from consecutive blocks of pulses, into a final fluctuation spec- 
trum. This can overcome the problem of low S/N observations. 
However, for pulsars which are affected by phenomena like drift 
mode changing, nulling or other effects, the results of the afore- 
mentioned spectral analysis do not show pure drift features. The 
process of averaging also results in a loss of temporal infor- 
mation, making the detection and characterisation of any time- 
related changes very difficult. In this section we propose an 
extension to the 2DFS technique which can resolve temporal 
changes of drift rates of subpulses. 

We call this method the Sliding Two-Dimensional 
Fluctuation Spectrum (S2DFS). The spectrum is obtained 
in the same way as in the 2DFS but instead of dividing a pulse 
stack, which consists of N pulses, into consecutive blocks of 
equal length, the 2DFS is applied to a block of M pulses, and 
a spectrum is calculated. The spectrum is collapsed horizon- 
tally producing a one-dimensional curve (see for example the 
left-hand side panel of the 2DFS plot in Fig.|2l right-hand plot). 
Subsequently, the DFT "window" is shifted by one pulse and 
the whole process is repeated. The application of the sliding 
DFT window to the pulse stack results in - M -i- 1 curves, 
which are stacked together to produce a so-called "map" of 
the collapsed fluctuation spectra (see Fig. |2] left-hand plot). 
Changing the position of the DFT window and calculating the 
2DFS allows one to resolve any changes of drift over relatively 
short timescales. The choice of the length of the DFT window 
is crucial for the resolution and sensitivity of the maps. If the 
size of sliding window, and hence the length of the Fourier 
transform, is too small, the resulting details will not be resolved 
due to the insufficient resolution in P^jP^ space, but if the size 
of the window is too large then the method will not be sensitive 
enough, and any short-lasting events visible in the pulse stack 
will not be sufficiently prominent in the maps due to the reduced 
S/N per spectral bin. In our work we have decided to use a 
DFT window size of 256 pulses and as will be shown later it 
has proven to be a good compromise between resolution and 
sensitivity. We note that in order to generate the aforementioned 
maps of the collapsed fluctuation spectra one can also use the 
LRFS method. We have decided to use 2DFS-based procedure, 
because collapsing the fluctuation spectra vertically can be 
used to obtain the information regarding changes of P2 values. 
However, the maps of vertically collapsed fluctuation spectra, 
holding information about the f 2 temporal changes, usually 
leads to results comparable to those from the 2DFS method 
alone. 



3. Simulated data 

To test the new technique of subpulse modulation analysis we 
have generated a number of artificial pulse sequences. We used 
an emission model where the emission is generated by a ro- 
tating "carousel" of discharges, or so-called "sparks", which 
circulate around the magnetic axis due to an E x B drift 
jRuderman & Sutherlandl[l975h . We know there are many prob- 
lems with the carousel model, but it proves to be sufficient for 
generating pulse stacks for testing our new algorithm, and we are 
not at this stage, concerned about the origin of the mechanism 
that produces the drifting subpulses. We have chosen the set of 
elementary model parameters used to generate the pulse stacks 
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Coherent drifting 




500 1000 1500 Fluctuation frequency (P„/P2) 

Block number 



Fig. 2. The analysis results shown for a simulated pulse stack assuming the coherent drifting scenario. Left: The S2DFS. The vertical 
axis is given in P0/P3 (cpp) and corresponds to the vertical separation of drift bands. The horizontal axis is given in blocks, where 
each block is a curve obtained from the vertically collapsed 2DFS where the DFT window of 256 pulses "slides" along the pulse 
stack. Subsequently, all curves are stacked together producing maps of the collapsed fluctuation spectra. Right: The fluctuation 
analysis results. The top panel shows the integrated pulse profile (solid line), the longitude-resolved modulation index, LRMI (solid 
line with error bars). The middle panel shows the LRFS where the ordinate of the resulting spectrum is given in cycles per period 
(cpp) which corresponds to P0/P3 (Po is the pulsar period and f 3 denotes the vertical separation of drift bands) and the abscissa 
denote pulse phase given in degrees. The bottom panel represents the 2DFS where the ordinate is the same as the LRFS but the 
abscissa is in cpp (^0/^2) (Pi denotes the horizontal separation of drift bands). The greyscale intensity of the 2DFS corresponds to 
the spectral power. The presence of a significant spectral feature with a negative or positive value of f 2 indicates that the subpulses 
appear with a preferred direction of drifting subpulses with periodicity of f pulsar periods. The side panel corresponds to the 
horizontally (left-hand panel) and vertically (bottom panel) integrated spectrum between the dashed lines. 



to be similar to the pulsar PSR B0809-(-74. These parameters are 
kept identical in all the pulse stacks in order to make comparison 
between the data sets easier Specifically, we used the following 
pulsar geometry: the angle between rotation and magnetic axes 
of the pulsar, a = 9 deg and the angle between rotation axis 
and the line of sight, ^ - 13.5 deg. The pulse stacks consist of 
a total of 2048 pulses with a period of Pq = Is. The period 
value is not used in our calculations and has no importance for 
the analysis presented in this paper We set the number of sparks 
N = 5 and we have chosen a drift rate value f 3 = 10 Po- Each 
subpulse has a full width at half-maximum (FWHM) of W50 = 5 
deg. The pulse period is divided into one thousand longitude bins 
and the average pulse profile has a FWHM of wso -12 deg and 
S/N= 100. 

The pulse sequences have been arranged into a few types of 
emission scenarios to investigate the effects of nulling, changing 
drift rates and mode changing on the S2DFS. In addition, to each 
emission scenario we introduced several effects which alter the 
fluctuation analysis results. The first effect was pulse-to-pulse in- 
tensity modulation with a normally distributed range of subpulse 
amplitudes. This makes the modulation pattern less pure, hence 



it broadens the feature in the spectra, making it more difficult 
to detect. The second effect was a steady emission component 
with a FWHM of W50 = 12 deg and an amplitude of 10 % of 
the average pulse profile. In some cases we also combined both 
the steady emission component and pulse-to-pulse modulation 
thogether 

3.1. Coherent drifting 

The first emission scenario is a simulation of coherently drifting 
subpulses with a chosen drift value of P3 - 10 Po- The 2DFS of 
this scenario (see Fig. |2l right-hand plot) showed a very strong 
and narrow (< 0.05 cpp) feature, at 0.1 cpp, corresponding to a 
P3 = 10 Po with P2 - 20 cpp (corresponding to 360 deg / 20 
cpp =18 deg pulse longitude). As expected, as there were no 
distorting effects, the S2DFS plot shows no change of drift rate 
over time in any of the simulated sequences. 

The introduction of pulse-to-pulse modulation, a steady 
emission component or both effects combined did not strongly 
affect the spectral results (not shown). However, after introduc- 
ing pulse-to-pulse modulation we noted the presence of a very 
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Drift change with nulling 
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Fig. 3. The analysis results shown for a simulated pulse stack assuming a drift rate change in combination with nulling. For the 
explanation of the plots we refer to Fig. |2] and the text. 



Table 1. The details of the analysed pulsars. The last two columns hold results presented in the work of I Weltevrede et akl (l2006l 
I2007h . 



Name 


Date V 


Po 


^ pulses 


Pi 


Pi 


PSR 


[MJD/dd.mm.yy] [MHZ] 


[s] 




[PiA 


[deg] 



B0031-07 


51780/24.08.2000 


328 


0.9430 


7666 


6.7 + 0.1 


-21 9+"^ 




53239/22.08.2004 


1380 




1204 


8.3 + 0.3 


-40+^ 

^"-50 


B1819-22 


53182/26.06.2004 


328 


1.8743 


1095 


16.9 + 0.6 


— 




53188/02.07.2004 


1380 




1096 


9.8 + 0.6 


-9 r" ' 


B 1944+ 17 


52706/08.03.2003 


328 


0.4406 


3990 


13.5 + 0.9 


-"-'-10 




51977/09.03.2001 


1380 




2056 


13 + 0.5 


-11 7+" 



weak, low-frequency feature around P0/P3 and P0/P2 = in 
the 2DFS plots of the analysed pulse stacks. We also note that 
despite the relatively low S/N all of the above spectral results 
showed second and third harmonics of the main drift feature in 
this scenario (see Fig.|2]i. We interpret this as being due to the co- 
herence and non-si nusoidal shape of the drift bands as discussed 
in E dwards & Stapp ers (2002). 

In the next scenario, we introduced nulling to the pulse stack. 
The null lengths and positions were randomly generated with a 
minimum nulling length of 5 pulses and total nulling fraction, 
NF = 22.9%. As in the coherent drifting scenario, the 2DFS 
showed a strong and narrow feature with P^ = 10 Po- We also 
note the presence of a low-frequency feature in the 2DFS plots 
(not shown), which we identify as being due to the presence of 
nulls. This feature is similar to the feature corresponding to the 
presence of pulse-to-pulse intensity modulation. We note that 
nulling, due to its timescales and occurrence rates, can be seen 



as an extreme form of intensity modulation, hence it is seen 
as a feature close to the P0/P3 - border. Also, due to the 
nulling only the second harmonic was visible in both the 2DFS 
and S2DFS plots in this scenario due to the weaker coherence of 
the modulation pattern. 



3.2. Drift rate change 

The next two scenarios included a constant change of drift rate, 
where the initial value of P3 =8.5 Pq changes to a value of 
P3 = 4.5 Po after 2048 periods, and in the case of the second 
scenario nulling was also introduced with a nulling fraction of 
NF = 22.9%. In both cases the 2DFS plots showed eight closely 
located drift features spanning 0.1 cpp, along the P0/P3 axis 
and in the second scenario a low-frequency feature as in the co- 
herent drift with nulling scenario is seen (see Fig. |3] right-hand 
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Null— induced mode change 




500 1000 1500 Fluctuation frequency (P„/P2) 

Block number 

Fig. 4. The analysis results shown for a simulated pulse stack assuming null-induced mode changes. For the explanation of the plots 
we refer to Fig. |2]and the text. 



plot). These unusual looking spectra are caused by the constant 
change of P3 and a property by which the 2DFS is calculated. 
This property is related to the process of calculating and averag- 
ing multiple spectra to obtain the final spectrum. When the pulse 
stack is divided into blocks of pulses and the DFT is applied to 
each consecutive block, each block's spectrum shows spectral 
features offset in the vertical direction from the spectrum of the 
first block. Those features correspond to the different values of 
the changing P3. The resulting block spectra were averaged to- 
gether to form the final 2DFS plot. What determines the number 
of spectral features in the 2DFS is the number of pulses in the 
pulse stack and the size of the block of pulses used in the DFT. 



In the case of the low S/N data it is very difficult to inspect 
the pulse stack by eye, making the spectral analysis the only re- 
liable method to determine the initial drift rate or the speed of 
its change. However, from the 2DFS alone one is not able to de- 
termine such drift rate variations. The S2DFS technique allows 
one to circumvent this problem as the process of averaging is re- 
placed by a moving DFT window. In the left-hand plot of Fig. [3] 
an example S2DFS plot is presented where one can easily see 
the "tracks" corresponding to the changing value of the main 
drift feature and its second harmonic. By analysing the S2DFS 
plot it is quite straightforward to determine the initial value of 
f 3 and its rate of change. Also inspection of the S2DFS plot can 
reveal the presence of nulling. This can be seen in the S2DFS 
plot (Fig. 12 left-hand plot) between blocks 100 to 200 or 1400 
to 1500 at 0.02 cpp. 



3.3. Null-induced mode change 

For the final simulated pulse sequence we investigated the null- 
induced mode change scenario. The generated pulse stack had 
two drift modes, mode "A" with P3 = 10 Pq s, and mode "B", 
P3 = 6 P() s, where the numbers of pulses in each mode was ran- 
domly chosen. The mode change was induced after a sequence 
of null pulses. The null lengths were also randomly generated. 
The minimal nulling length which triggered the change of the 
drift modes was set to eight periods, the overall nulling frac- 
tion was, NF = 32.4% and the calculated numbers of pulses in 
the drift modes "A" and "B" contained 797 and 587 pulses re- 
spectively. The return to the preceeding drift mode was executed 
abruptly without recovery time. 

Fig.|4](bottom right-hand plot) presents the 2DFS of this case 
where a pulse-to-pulse modulation and a steady emission com- 
ponent were also introduced to the pulse stack. In the spectrum 
one can easily see many features present between Po/Pj = 0.05 
and 0.35 cpp. In the side panel showing the horizontally col- 
lapsed spectrum, the peak at P0/P3 = 0.05 cpp is the most 
prominent and we identify it as caused by the nulls. The next 
two peaks at 0. 1 and 0. 166 cpp we identify them as the main drift 
features (P3 = 10 f and Pj, - 6 Po respectively), whereas the 
pair of peaks visible at P0/P3 - 0.2 and 0.33 cpp are their sec- 
ond harmonics. Inspection of the S2DFS plot (Fig. |4] left-hand 
plot) allows us to determine the "A" mode as being the most 
prominent. One can easily see the change from the "A"-mode 
to "B"-mode around block 1100. We also note that in first 700 
blocks both drift modes are not easily seen in the S2DFS plot. 
We interpret this as being due to the frequent changes between 
drift modes, which resulted in a decrease of the spectral power 
of features corresponding to the drift modes. The S2DFS plot 
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Fig. 5. The analysis results shown for PSR B003 1 -07. For the explanation of the plots we refer to Fig. |2]and the text. 



also explains the remaining features present around 0. 1 cpp (e.g. 
block 400) as being created during the abrupt switch between 
the drift modes. 



4. Observations 

The above results clearly present the S2DFS technique as a com- 
pliment of the 2DFS method for detecting and analysing the 
drifting subpulses phenomenon, particularly for sources with 
multiple or changing drift rates. To test the new method us- 
ing real observations, we have selected a number of archival 
data sets taken at the 94-m equivalent Westerbork Synthesis 
Radio Telescope (WSRT) in the Netherlands. The three selected 
pulsars, PSRs B0031-07, B1819-22 and B1944+17, were ob- 
served between 2000 and 2004 at two wavelengths of 21 and 92 
cm. Specifically, we have used the PSR B003 1 -07 data from the 
work of Smits et al. (2005), while the data of PSR B1819-22 
and PSR B 1944+ 17 were taken from IWeltevrede et all (l2006i 
l2007i) . We have selected these sources based on the quality of 
data and the well known properties of these pulsars. This will 
allow us to compare the results from the 2DFS analysis with the 
new S2DFS technique. For a detailed disc ussion of the obser ving 
system and data reducti on procedures see lVoute et al.l(l2002l) and 
IWeltevrede et al.l(l2006l) . while Table[T|summarises the details of 
the obse rving sessions and results from the 2DFS analysis pre- 
sented in lWeltevrede et al.1 (I2006il2007l) . 

4.1. PSR B0031-07 

This pulsar shows a broad drifting feature in its 2DFS and 
S2DFS at 92 cm (Fig.|5] left-hand and middle plot). In the case of 
the S2DFS plot we only show first 1500 blocks in order to keep 



the figure legible as the whole data set has 741 1 blocks. Three 
drift modes have been found for this pulsar by Hugueni n et al. 
(11970.) at 145 and 400 MHz. In the work of .Weltevrede et al. 
(I2007h most of the power in the 2DFS at 92 cm is due to the "B"- 
mode drift (P3 = 6 f 0) at 0.15 cpp. The slope of the drift bands 
chang es slightly from band to band (e.g. Vivekanand & Josh| 
1 1997b . causing the feature to extend vertically in the 2DFS plot. 
The results from the S2DFS analysis confirm the "B"-mode drift 
as the most prominent one. One can see that the region, in the 
S2DFS plot, corresponding to the "B"-mode drift seems to con- 
sist of three separate tracks. We interpret this as being caused by 
the changing slopes of the drift bands. 

The "A"-mode drift (P3 = 12 Pq) is visible as a downward 
extension of the main drift feature to 0.08 cpp as can be seen in 
the 2DFS greyscale plot (Fig. |5] bottom panel of middle plot). 
During the "A"-mode the pulsar signal becomes weaker which 
can be directly seen in the pulse stack and as a decrease of spec- 
tral power of the feature in the S2DFS plot. In the case of the 
S2DFS an exact location and duration of the track correspond- 
ing to this drift mode can be seen in the left-hand plot of Fig. |5] 
between blocks 750 and 1000 at 0.08 cpp. Two more regions 
(not shown), denoting the presence of the "A"-mode drift, were 
revealed after visual inspection of the whole S2DFS plot for 
this pulsar. We note the presence of a low-frequnecy feature, 
located close to P0/P3 = border. We identify both pulse-to- 
pulse modulations and nulling responsible for the presence of 
this feature as confirmed in the results from the analysis of sim- 
ulated data. However, during the presence of the "A"-mode drift 
the low-frequency feature disappears. We interpret this as fol- 
lows: the region in the pulse stack corresponding to the presence 
of the slow drift mode shows no nulling as well as a decrease 
of the modulation index reaching the value of 1.3, which in- 
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Fig. 6. The analysis results shown for PSR B 18 19-22. For the explanation of the plots we refer to Fig.|2]and the text. 



dicates low pulse-to-pulse modulation. The "C"-mode drift at 
0.25 cpp (P? ^ - 19+? deg and P3 = 4.1 + 0.1 Po) reported by 
IWeltevrede et~ai] (12007) as seen in the 2DFS is also visible, al- 
though weakly, in the S2DFS plot (not shown). In both methods 
the feature at 0.3 cpp is the second harmonic of the main drift 
feature. 

In the 21-cm data, IWeltevrede et al.1 (120061) report the "A"- 
mode dominated over the "B"-mode which is visible, although 
weakly, while the "C"-mode was not detected in the 2DFS. In 
the case of results using the S2DFS (not shown for the 21 cm 
data) only the "A"-mode is seen in the plot. The remaining "B"- 
and "C"-drift modes were not detected as the insufficient S/N 
makes it difficult to discriminate the features responsible for the 
drift from the noise in the S2DFS plot. The behaviour of this 
pulsar in the 92 cm and 21 cm observations i s consistent with th e 
multifrequency study of lSmits etaP (l2005h : ISmits et all (l2007h . 

4.2. PSR B1 81 9-22 

The results from the 21 cm data show a feature in the 2DFS 
plot (Fig. |6] bottom panel of the central plot) which is double- 
peaked and corresponds to the two separate drift modes. The 
slower drift mode, present at . 06 cpp and fast drift mode present 
at 0.16 cpp (I Weltevrede et all 12006'). In Fig. |6] (left-hand plot), 
the S2DFS plot clearly shows that the fast drift mode appears 
between blocks 300 and 600 at 0.16 cpp, corre sponding to P3 - 
8 ± 4 f and P2 - -50+[g° deg as reported bv IWeltevrede etalj 
(I2OO7I) (see Table [Til. We note that the visual inspection of the 
pulse stack revealed sparsely present fast mode drift bands in 
the whole pulse stack. In the case when only solitary drift bands 
are present one cannot see any contribution in the S2DFS, while 
just few consecutive drift bands can be resolved by the S2DFS 



as can be seen in Fig. |6] (left-hand plot, e.g. blocks 80 to 100). 
We also note that the drift mode change is preceded by a null, 
that can last from few down to a single pulse period and can be 
seen only after a visual inspection of the pulse stack. 

We note several regions, in the S2DFS plot, where one can- 
not see clear tracks denoting drift (Fig. |6] (left-hand plot, e.g. 
blocks 100 to 170), we interpret this as either the presence of 
nulling or rapidly changing drift modes causing smearing out in 
the DFT because a too long transform length was used. However, 
introducing shorter transform length decreased the spectral res- 
olution of the spectra and changed the S/N per spectral bin. This 
made the identification of features marking rapidly changing 
drift modes unattainable. 

The 2DFS and S2DFS of this pulsar at 92 cm very clearly 
show a drift feature at 0.06 cpp (Fig.|6] bottom panel of the right- 
hand plot) which corresp onds to P3 - 16.9 + 0.6 Pq and P2 = 
-17^g deg as reported bv IWeltevrede et aTl (l2007h . This feature 
corresponds to the feature denoting the slow drift mode present 
in the 21 cm data while the fast drift mode is not present in the 
92 cm data. There is a low-frequency excess present in the 2DFS 
plot (Fig.|6l bottom panel of the right-hand plot) and it is stronger 
than in the 21 cm data. We interpret this as being due to the 
frequent nulls appearing in the pulse stack instead of a fast drift 
mode pulses. 

4.3. PSR B 1944+ 17 

This pulsar shows drift features in its 2DFS at 92 and 21 cm 
(Fig- 111 middle and right-hand plot) and the drifting can clearly 
be seen by eye in pulse stacks at both observed wavelengths. 
The features in the 2DFS pl ot are broadened b ecause this pulsar 
shows drift-mode changes (iDeich et al.lll986l at both 430 and 
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Fig. 7. The analysis results shown for PSR B 1944+ 17. For the explanation of the plots we refer to Fig.|2]and the text. 



1420 MHz). The P3 = 13 Po "A"-mode drift is visible in the 
2DFS plot at 0.075 cpp at 92 and 21 cm. 

The results from the S2DFS analysis for both wavelengths 
showed the presence of tracks in the S2DFS plot revealing the 
presence of the aforementioned drift mode. The sparse occur- 
rence and very short duration times of the tracks (of the order of 
one hundred blocks) were caused by the frequent occurrence of 
nulls disrupting the drift bands. In the left-hand plot of Fig.|2]we 
show the S2DFS plot of 1300 out of 3735 blocks with the most 
prominent example of tracks denoting "A"-mode drift present at 
92 cm. 

One can easily see that between block 3100 and 3150 there 
is a region where power in the S2DFS plot is dominated by fluc- 
tuations across the whole Po/Pj, range. The visual inspection of 
the pulse stack revealed the presence of a group of pulses, which 
did not show any particular organised order, in between two long 
nulls disrupting the "A"-mode drift. We interpret this as follows: 
The DFT window sliding through the long null, enclosed the 
group of pulses, and thus produced fluctuations across the Pq/Pj 
range in the collapsed fluctuation spectra. One can also see that 
after block 3500 the low-frequency feature decreases its spectral 
power. We explain that as being caused by similar behaviour as 
in PSR B0031-07. 

We also note that after this null the position of the track in the 
S2DFS plot denoting the "A"-mode drift (Fig. |2l left-hand plot, 
block 3 150 to 3450) changes its position with respect to the track 
before the null. The 2DFS analysis of the pulse stacks containing 
the drift bands from the aforementioned region revealed that the 
drift rate value changes from P3 = 13.5 + 0.9 Po to P3 = 14.3 + 
0.6 P(). Without the use of the S2DFS it would not be possible to 
see that these changes are discrete changes in dirft rate and not 



smooth. It also shows that PSR B 1944+ 17 belongs to the class 
of null-induced drift mode changing pulsars. 

The P3 = 6.4 Po "B"-mode drift (0.16 cpp) did not appear 
in the 2DFS plot as a peak, although the centroid of the power 
was offset from the vertical axis by up to at least 0.2 cpp in the 
vertical direction at the wavelength of 21 cm (see Fig. |7] right- 
hand plot) indicating its presence is either weak or only very 
sporadic. This drift mode was only visible in the S2DFS plot in 
the 2 1 cm observation. We note the low-frequency excess in both 
observations due to the presence of nulling. 



5. Summary and conclusions 

We have presented a new technique called the Sliding Two- 
Dimensional Fluctuation Spectrum (S2DFS), which is a good 
supplementary tool to the 2DFS and is capable of discovering 
and characterising the temporal changes of drifting subpulses 
from radio pulsars caused by e.g. null-induced mode or drift 
rate changes, something which cannot be done using the 2DFS 
method alone. We note that for the S2DFS method the S/N of 
the analysed data and length of the DFT window play a crucial 
role in obtaining sufficient quality of the maps of the collapsed 
fluctuation spectra. We have chosen the size of the DFT window 
to be 256 pulses which proved to be the best compromise be- 
tween the resolution of the spectra (and hence the S2DFS plot) 
and sensitivity of the method to the short-lasting events like e.g. 
mode changes or nulls. 

In order to test the S2DFS method we have generated a num- 
ber of artificial pulse stacks. The emission models comprised of 
five scenarios: coherently drifting subpulses, coherently drifting 
subpulses with nulling, drift rate change, drift rate change with 
nulling and null-induced mode change. The set of pulse stacks 
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generated for testing purposes comprised of the aforementioned 
different emission scenarios with additional effects like steady 
emission component or a pulse-to-pulse modulation introduced 
to the pulse stacks. These pulsar-intrinsic emission effects are 
known to disrupt drift bands and alter the results of spectral anal- 
ysis. The maps of the collapsed fluctuation spectra allowed us 
to recover occurrence times of mode changes in the case of the 
null-induced mode change scenario and we were also able to 
resolve which one of the modes was more prominent. The sim- 
ulated drift rate changes were easily revealed, something that 
could not have been done by the 2DFS only. 

We have also used data sets from archival WSRT observa- 
tions at wavelengths of 92 and 21 cm. We have analysed three 
pulsars, PSRs B0031-07, B1819-22 and B1944+17, which 
were selected based on the quaUty of the data and their known 
drift properties. All three sources are known to exhibit mode 
changes which is most prominent in the case of PSR B 1 8 1 9-22 
at the wavelength of 21 cm and could be easily seen in the 
S2DFS plot. The usefulness of the new method is especially ap- 
parent in the case of short periods of different drift modes, where 
it is easier to find the moment of change, relative strength or co- 
herence of the drift mode without having to visually inspect the 
pulse stack. Moreover, inspection of the pulse stack does not di- 
rectly give the drift rate values. By applying this method to a 
large sample of pulsar such analysis would help to understand 
when and why drifting is non-coherent in different sources. We 
also note the potential of this method when applied to milisecond 
radio pulsars. If drift rate changes are discovered in recycled pul- 
sars it would give significant insight into the pulse profile evo- 
lution of miUisecond pulsars. This knowledge could be used for 
the construction of more precise timing models and hence, im- 
proving timing precision. 
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